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Abstract 

We consider a weak scale super symmetric seesaw model where the 
Higgsino is the next-to-lightest supersymmetric particle and the right- 
handed sneutrino is the dark matter candidate. It is shown that, in 
this model, inverse decays, which had been previously neglected, may 
suppress the sneutrino relic density by several orders of magnitude. 
After including such processes and numerically solving the appropriate 
Boltzmann equation, we study the dependence of the relic density 
on the /i parameter, the sneutrino mass, and the neutrino Yukawa 
coupling. We find that, even though much smaller than in earlier 
calculations, the sneutrino relic density is still larger than the observed 
dark matter density. 



1 Introduction 



The supersymmetric version of the seesaw mechanism is an attractive candi- 
date for physics beyond the Standard Model. On the one hand, it includes the 
seesaw mechanism, which postulates the existence of right-handed neutrino 
fields and has become the most popular framework to account for neutrino 
masses. The seesaw is able to accommodate the experimental data on neu- 
trino masses and mixing pLj, explaining naturally the small neutrino mass 
scale. On the other hand, it embraces low energy supersymmetry, with its 
rich phenomenology and its well known virtues. In fact, the minimal su- 
persymmetric Standard Model solves the hierarchy problem, achieves the 
unification of the gauge couplings, and contains a dark matter candidate: 
the lightest supersymmetric particle. 

The lightest sneutrino is a new dark matter candidate present in the 
supersymmetric seesaw. Being a mixture of left-handed and right-handed 
sneutrino fields, the lightest sneutrino will have different properties depend- 
ing on its composition in terms of interactions eigenstates. In general, three 
different kind of sneutrinos can be envisioned: a dominantly left-handed one, 
a mixed sneutrino, or a dominantly right-handed one. A dominantly left- 
handed sneutrino is not a good dark matter candidate. They are ruled out 
by experimental searches [2] and tend to have a too small relic density [3]. A 
mixed sneutrino can be compatible with the observed dark matter density as 
well as with present bounds from direct searches [31[5]. The required mixing 
is obtained at the expense of a large neutrino trilinear coupling, which is not 
allowed in typical models of supersymmetry breaking. A dominantly right- 
handed sneutrino is the final possibility, the one we will be concerned with 
throughout this paper. A right-handed sneutrino, being essentially sterile, 
interacts with other particles mainly through the neutrino Yukawa coupling. 
Could such a sterile sneutrino account for the observed dark matter density? 

Gopalakrishna, Gouvea, and Porod, in [6], studied that possibility within 
the same scenario we are considering here. They showed that self-annihilations 
of right-handed sneutrinos as well as co-annihilations with other particles are 
too weak to keep the sneutrinos in equilibrium with the thermal plasma in 
the early Universe. They also found that the production of sneutrinos in 
the decay of other supersymmetric particles gives a too large contribution 
to the relic density. They concluded, therefore, that in the standard cosmo- 
logical model right-handed sneutrinos cannot explain the dark matter of the 
Universe. 
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Even though generally valid, that conclusion is not guaranteed if the 
mass difference between the Higgsino and the sneutrino is small. In that 
case, inverse decays, such a.s N + L ^ H, contribute to the annihilation of 
sneutrinos and therefore to the reduction of the sneutrino relic density. Such 
possibility was not taken into account in |6|. In this paper, we will focus 
on models with a Higgsino NLSP and show that inverse processes cannot 
be neglected, for they suppress the sneutrino relic density by several orders 
of magnitude. Then, we will reexamine whether the sterile sneutrino can 
explain the dark matter of the Universe in the standard cosmological model. 

In the next section we briefly review the supersymmetric seesaw model 
and show that sterile sneutrinos arise naturally in common scenarios of su- 
persymmetry breaking. Then, in section [31 we will include inverse decays 
into the Boltzmann equation that determines the sneutrino abundance. It 
is then shown that inverse decays are indeed relevant; they cause a signifi- 
cant reduction of the relic density. In section HI we study the relic density 
as a function of the neutrino Yukawa coupling, the sneutrino mass, and the 
Higgsino-sneutrino mass difference. There, we will obtain our main result: 
the suppression effect of inverse decays, though important, is not enough to 
bring the sneutrino relic density down within the observed range. In the final 
section we will review our study and present our conclusions. 

2 The model 

We work within the supersymmetric version of the seesaw mechanism, where 
the field content of the MSSM is supplemented with a right-handed neutrino 
superfield N per generation. The superpotential then reads 

W = Wmssm + ^M'^'n'n' + YJ'H^L'n' (1) 

where, as usual, we have assumed R-parity conservation and renormalizabil- 
ity. Mn is the Majorana mass matrix of right-handed neutrinos and Y^y is 
the matrix of neutrino Yukawa couplings. Without loss of generality Mn 
can be chosen to be real and diagonal. Yjy is in general complex but we will 
assume, for simplicity, that it is real. Mn and Yi, are new free parameters 
of the model; they are to be determined or constrained from experimental 
data. 



2 



After electroweak symmetry breaking, the above superpotential generates 
the following neutrino mass terms 



i:-umass = -VuY.uN - ^M^NN + h.c. (2) 

If M^v ^ v^Y,^, the light neutrino mass matrix, rriy, is then given by the 
seesaw formula 

= -mnM^^m^, (3) 

with rriD = VuY^ being the Dirac mass. Since mu is partially known from 
neutrino oscillation data, equation ([3]) is actually a constraint on the possible 
values of Y^ and M^y. It is a weak constraint though; and it allows to 
vary over many different scales. In this paper we consider what is usually 
known as a seesaw mechanism at the electroweak scale. That is, we assume 
that Mtv ~ 100 GeV. Thus, since the neutrino mass scale is around nii, ~ 0.1 
eV, the typical neutrino Yukawa coupling is 

Yu ~ 10-6 , (4) 

or around the same order of magnitude as the electron Yukawa coupling. 
Notice that this value of Y^ is a consequence of the seesaw mechanism at the 
electroweak scale. In other frameworks, such as Dirac neutrinos or seesaw 
at much higher energies, Y^ takes different values. We will not consider such 
possibilities here. 

The new soft-breaking terms of the supersymmetric seesaw model are 
given by 

Csoft = -{miy'N*^Ni + [{mlY'NiNi - A'^KUNi + h.c] . (5) 

They include sneutrino mass terms as well a trilinear interaction term. For 
simplicity, we will assume that m%, and Ay are real. 

To study the sneutrino mass terms resulting from ([1]) and ([5]) it is conve- 
nient to suppress the generation structure; that is, to work with one fermion 
generation only. It is also useful to introduce the real fields i>i, z>2, A''i and 
N2 according to the relations 

i>L = (i>i + iz>2) ,Nr = ^ (iVi + iN^ . (6) 
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Indeed, in the basis {ui, Ni, 1/2, N2) the sneutrino mass matrix takes a block 
diagonal form 



I m\^ + mj:,Mi^ 



- mnM, 



N 



\ ^RL ~ ^dMn f^RR + 



(7) 

where rriLL = ^| + + 0.5m| cos 2/5, m\p^ = M% + m\ + m|,, and m\j^ = 
—fiVdY^ + f uv4^. This matrix can be diagonalized by a unitary rotation with 
a mixing angle given by 

tan26'i''2 = — — ^ — tt^ (8) 

' mlL-i^RRTml)' 

where the top sign corresponds to 61 -to the mixing between ui and A^i- 
whereas the bottom sign corresponds to ^2- 

Since A^,> is independent of gaugino masses, there is a region in the 
super symmetric parameter space where the lightest sneutrino, obtained from 
(I7j), is the lightest supersymmetric particle (LSP) and consequently the dark 
matter candidate. That is the only region we will consider in this paper. 

The lightest sneutrino is a mixture of left-handed and right-handed sneu- 
trino fields. Depending on its gauge composition, three kinds of sneutrinos 
can be distinguished: a dominantly left-handed sneutrino, a mixed sneutrino, 
and a dominantly right-handed sneutrino. A dominantly left-handed sneu- 
trino is not a good dark matter candidate for it is already ruled out by direct 
dark matter searches. These sneutrinos also have large interactions cross sec- 
tions and tend to annihilate efficiently in the early universe, typically yielding 
a too small relic density. A mixed sneutrino may be a good dark matter can- 
didate. By adjusting the sneutrino mixing angle, one can simultaneously 
suppress its annihilation cross section, so as to obtain the right relic density, 
and the sneutrino-nucleon cross section, so as to evade present constraints 
from direct searches. A detailed study of models with mixed sneutrino dark 
matter was presented recently in [5J. A major drawback of these models is 
that the required mixing may be incompatible with certain scenarios of su- 
persymmetry breaking, such as gravity mediation. The third possibility, the 
one we consider, is a lightest sneutrino which is predominantly right-handed. 
That is, a sterile sneutrino. 
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A sterile sneutrino is actually unavoidable in supersymmetry breaking 
scenarios where the trilinear couplings are proportional to the corresponding 
Yukawa matrices, such as the constrained Minimal Supersymmetric Standard 
Model (CMSSM)[I]. In these models 

= a^X^irLsoft (9) 

where nisoft ~ 100 GeV is a typical supersymmetry breaking mass and is 
an order one parameter. Because Y^, is small, A^, is much smaller than the 
electroweak scale, 

A^-lOOkeV. (10) 

Hence, from equation ([HI), the mixing angle between z/j and iVj is also very 
small 

sinei~10-^ (11) 

Thus, we see how in these models the small Yy translates into a small trilinear 
coupling Ay that in turn leads to a small mixing angle -to a sterile sneutrino. 
Sterile sneutrinos are also expected in other supersymmetry breaking mech- 
anisms that yield a small A^ at the electroweak scale. 

Since the mixing angle is small, we can extract the sterile neutrino mass 
directly from ([7]). It is given by 

^\ = ^RR-^B-Ml + m\-m\ (12) 

where we have neglected the Dirac mass term in the last expression, 

is thus expected to be at the electroweak scale. In the following, we will 

consider m^y = "^l5p as a free parameter of the model. 

To summarize, the models we study consist of the MSSM plus an elec- 
troweak scale seesaw mechanism that accounts for neutrino masses. Such 
models include a new dark matter candidate: the lightest sneutrino. In com- 
mon scenarios of supersymmetry breaking, the lightest sneutrino, which we 
assume to be the dark matter candidate, turns out to be a dominantly right 
handed sneutrino, or a sterile sneutrino. In the following, we will examine 
whether such a sterile sneutrino may account for the dark matter of the 
Universe. 

3 The N relic density 

To determine whether the sterile sneutrino can explain the dark matter of 
the universe we must compute its relic density fi^/i^ and compare it with 



5 



the observed value VL^Mh'^ = 0.11[H]. This question was aheady addressed in 
[B]. They showed that, due to their weak interactions, sneutrinos are unable 
to reach thermal equilibrium in the early Universe. In fact, both the self- 
annihilation and the co-annihilation cross section are very suppressed. They 
also noticed that sneutrinos could be produced in the decays of other super- 
symmetric particles and found that such decay contributions lead to a relic 
density several orders of magnitude larger than observed. Thus, they con- 
cluded, sterile sneutrinos can only be non-thermal dark matter candidates. 

That conclusion was drawn, however, without taking into account inverse 
decay processes. We now show that if the Higgsino-sneutrino mass difference 
is smal0, inverse decays may suppress the sneutrino relic density by several 
orders of magnitude. To isolate this effect, only models with a Higgsino 
NLSP are considered in the following. We then reexamine the possibility 
of having a sterile sneutrino as a thermal dark matter candidate within the 
standard cosmological model. 

In the early Universe, sterile sneutrinos are mainly created through the 
decay H —>■ N + L, where H is the Higgsino and L is the lepton doublet. 
Alternatively, using the mass-eigenstate language, one may say that sneu- 
trinos are created in the decay of neutralinos ^ N + u) and charginos 
(x^ — >■ £^ + A^). These decays are all controlled by the neutrino Yukawa 
coupling Y^. Other decays, such as £ ^ Nff via W^, also occur but the 
Higgsino channel dominates. Regarding annihilation processes, the most im- 
portant one is the inverse decay N + L ^ H. In fact, the sneutrino-sneutrino 
annihilation cross section is so small that such process never reaches equi- 
librium. And a similar result holds for the sneutrino coannihilation cross 
section. We can therefore safely neglect annihilations and coannihilations in 
the following. Only decays and inverse decays contribute to the sneutrino 
relic density. 

The Boltzmann equation for the sneutrino distribution function then 
reads: 

^-ET^TE^ = TTET- / /o^AHoF /o„A3oz7 \^H^Ln\ (1'^) 



dt E dE 2Ej^ J {2'Kf2EL{2Tif2EH 



(2nr5\ps~PL-PN)[fH-fLl 



N\ 



where H is the Hubble parameter and ffj, fi respectively denote the H and 
L distribution functions. Other dark matter candidates, including the neu- 



^If it is large the results in [6] would follow. 
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tralino, have large elastic scatterings cross sections with the thermal plasma 
that keep them in kinetic equilibrium during the freeze out process. Their 
distribution functions are then proportional to those in chemical equilibrium 
and the Boltzmann equation can be written as an equation for the number 
density instead of the distribution function [9]. For sterile sneutrinos, on 
the contrary, the elastic scattering is a slow process -being suppressed by 
the Yukawa coupling- and kinetic equilibrium is not guaranteed. Hence, we 
cannot write (fT3l) as an equation for the sneutrino number density and 
must instead solve it for 

If the condition /jy ^ 1 were satisfied, inverse processes could be ne- 
glected and a simple equation relating the sneutrino number density to the 
Higgsino number density could be obtained. That is the case, for instance, in 
super symmetric scenarios with Dirac mass terms only [7]. In such models, the 
neutrino Yukawa coupling is very small, Y,y ~ 10~^^, and sneutrinos never 
reach chemical equilibrium. But for the range of parameters we consider, 

~ 10~^, the condition -C 1 is not satisfied. 

Since equation f[T^ depends also on the Higgsino distribution function, 
one may think that it is necessary to write the Boltzmann equation for 
and then solve the resulting system for /jy and fg. Not so. Higgsinos, due to 
their gauge interactions, are kept in thermal equilibrium -by self-annihilation 
processes- until low temperatures, when they decay into N+L through the Y^, 
suppressed interaction. It is thus useful to define a freeze-out temperature, 
Tf o,, as the temperature at which these two reaction rates become equal. 
That is, 

nH{crHHv)\T,.^, =TiH^N + , (14) 

where ng is the Higgsino number density and {crfjfjv) is the thermal av- 
erage of the Higgsino-Higgsino annihilation rate into light particles. Ty o. 
marks the boundary between two different regimes. For T > Tf o. Higgsinos 
are in equilibrium and annihilate efficiently. The Higgsinos produced in the 
inverse decays, in particular, easily annihilate with thermal Higgsinos into 
light particles. The inverse process is thus effective. In contrast, for T < Tj.o. 
Higgsinos mostly decay into the LSP and inverse decays cannot deplete the 
sneutrino abundance. The final state Higgsinos simply decay back into sneu- 
trinos: N + L-^H-^N + L. Below Ty.o., therefore, the total number of 
sneutrinos plus Higgsinos remains constant. Thus, we only need to integrate 
equation (ITSl) until Tf ^., a region in which Higgsinos are in equilibrium. 
Assuming a Maxwell-Boltzmann distribution, f{E) oc exp(— £'/T), for 
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Higgsinos and leptons and neglecting lepton masses, the integrals in (IT^ can 
be evaluated analytically to find 



dt 
where 



dt E dE 167rE^|p^| ^ ^ J V ; 



2 2 
m - — m- 



E^ = ^i-^{E^±\V^\). (16) 

In the following we will solve equation f|T5]) to obtain the sneutrino abun- 
dance, Yfj = fij^/s, and the sneutrino relic density, Qj^h'^. The sneutrino 
abundance today will be given by 

>^7vIto=>^^|t,.„. +r^|T,..., (17) 

where the second term takes into account that the Higgsinos present at freeze- 
out will decay into sneutrinos. The sneutrino relic density today is then 
obtained as 

a,h^ = 2.8 X io'"r^^. (18) 

The only parameters that enter directly in the computation of the sneutrino 
relic density are the Yukawa coupling, the sneutrino mass, and the Higgsino 
mass, which we take to be given by the fi parameter -fUj^ = All other 
supersymmetric particles besides N and H are assumed to be heavier, with 
msusy ~ 1 TeV. To determine the freeze-out temperature, equation f|T^ . we 
also need to know the Higgsino annihilation rate into Standard Model par- 
ticles. We use the DarkSUSY package [lU] to extract that value. Regarding 
the initial conditions, we assume that at high temperatures {T ^ rrifj) the 
sneutrino distribution function is negligible ff^r^O. Finally, we assume that 
the early Universe is described by the standard cosmological model. 

Once decays and inverse decays are included in the N Boltzmann equa- 
tion, two questions naturally come to mind. First, for what values of Y^, are 
inverse decays relevant? Second, can decays and inverse decays bring the 
sneutrinos into equilibrium? To answer these questions we show in figure 
[T]the sneutrino abundance as a function of the temperature for rriff = 100 
GeV, = 120 GeV, and different values of Fj.. Notice that for = 10"^ 
inverse processes are negligible and the sneutrino abundance simply grows 
with temperature. In that region, for Y^, < 10~^, the sneutrino relic density 
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Figure 1 : The role of the neutrino Yukawa couphng on the sterile sneutrino abun- 
dance. The figure shows y as a function of the temperature for different values of 
Yiy. The sneutrino mass is lOOGeV while n = 120GeV. 



is proportional to Y^^. From the figure we see that for Yi, = 10"^ the in- 
verse process leads to a reduction of the sneutrino abundance around T = 20 
GeV. The Yukawa interaction is not yet strong enough to bring the sneu- 
trinos into equilibrium. For Y^, = 10~^ sneutrinos do reach equilibrium and 
then decouple at lower temperatures. For even larger Yukawa couplings, 

= 10~^, 10~^, equilibrium is also reached but the decoupling occurs at 
higher temperatures. In that region, the relic density also increases with the 
Yukawas. Thus, for Y^ ~ 10~^ inverse decays not only are relevant, they are 
strong enough to thermalize the sneutrinos. 

Figure [2] directly compares the resulting sneutrino abundance with and 
without including the inverse process. The full line corresponds to the correct 
result, taking into account the direct and the inverse process. The dashed 
line, instead, shows the result for the direct process only, that is the sneutrino 
abundance according to [6]. The sneutrino mass was taken to be lOOGeV and 
Yi, was set to 10~^. The Higgsino mass is different in each panel and includes 
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Figure 2: The effect of the inverse process on the sneutrino relic density. The 
panels show the resulting sneutrino abundance Y = n/s as a function of the 
temperature for rrifj = lOOGeV and different values of fi. The full line is the result 
obtained including the inverse process whereas the dashed line is the result without 
including them. The dash-dotted line shows the sneutrino equilibrium abundance. 

values leading to strong and mild degeneracy as well as no-degeneracy at 
all between the sneutrino and the Higgsino. Notice that the correct final 
abundance, and consequently the resulting relic density, is always several 
orders of magnitude below the value predicted in p]. Even for the case of a 
large mass difference, we find a suppression of 3 orders of magnitude in the 
relic density. And as the mass difference shrinks the suppression becomes 
larger, reaching about 6 orders of magnitude for ^ = 150 and about 7 orders 
of magnitude for fi = 120GeV. We thus see that over the whole parameter 
space the inverse process has a large suppression effect on the sneutrino relic 
density. 
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Figure 3: The sneutrino relic density as a function of the neutrino Yukawa coupHng 
for different values of m^, and Am = 20GeV. 



4 Results 

So far we have found that the inverse decay process N + L ^ H leads to a 
suppression of the sneutrino rehc density. It remains to be seen whether such 
suppression is strong enough to bring the rehc density down to the observed 
value. That is, we will now study the dependence of the relic density with 
the sneutrino mass, the Higgsino-sneutrino mass difference, and the neutrino 
Yukawa coupling to find the region of the parameter space that satisfies the 
condition Qf^h"^ = VL^Mh?- 

Figure [3] shows the sneutrino relic density as a function of the neutrino 
Yukawa coupling and different values of the sneutrino mass. The Higgsino- 
sneutrino mass difference (Am = nifj — mf;) was set to 20 GeV. Larger values 
would only increase the relic density -see figure O Notice that, for a given 
sneutrino mass, the relic density initially decreases rather steeply reaching 
a minimum value at Fj, < 10^^ and then increases again. From the figure 
we also observe that the smallest value of the relic density is obtained for 
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Y„ = 1 .Oe-6 



"'too 150 200 250 300 350 

nYsp(GeV) 

Figure 4: The sneutrino relic density as a function of the sneutrino mass for 
yjy = 10~^ and different values of Am/m^y. As expected the smaller the mass 
difference the smaller the relic density is. 



TTifj = 400 GeV, that is, when the percentage mass difference is smaller. In 
any case, the relic density is always larger than 1, too large to be compatible 
with the observations. 

This result is confirmed in figure H] when we display the relic density 
as a function of the sneutrino mass for Yf^ = 10~^ and different values of 
Am/m. In agreement with the previous figure, we see that the smaller 
the percentage mass difference, the smaller the relic density is. Yet, Qfjh'^ is 
always larger than 1. We have verified that this conclusion is robust. Neither 
larger sneutrino masses nor different Yukawa couplings lead to the correct 
value of the relic density. 

5 Conclusions 

We studied the possibility of explaining the dark matter with a sterile sneu- 
trino in a supersymmetric model consisting of the MSSM supplemented with 
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a seesaw mechanism at the weak scale. We showed that if the Higgsino 
is the NLSP inverse decays play a crucial role in the computation of the 
sneutrino relic density, suppressing it for several orders of magnitude. We 
wrote down and numerically solved the correct Boltzmann equation that 
determines the sneutrino abundance and studied the resulting relic density 
as a function of the sneutrino mass, the neutrino Yukawa coupling and the 
Higgsino-sneutrino mass difference. We found that the sterile sneutrino relic 
density, even though much smaller than previously believed, is still larger 
than the observed dark matter density. In this scenario, therefore, the sterile 
sneutrino is not a thermal dark matter candidate. 
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